Hamiltonian Analysis of Lagrange Multiplier 

Modified Gravity 



Josef Kluson 



Department of Theoretical Physics and Astrophysics 
Faculty of Science, Masaryk University 
Kotldfskd 2, 611 37, Brno 
Czech Republic 



Abstract 

We develop Hamiltonian formalism for Lagrange Multiplier Modified Grav- 
ity. We further calculate the Poisson brackets between constraints and we 
show that they coincide with the algebra of constraints in Hamiltonian for- 
mulation of General Relativity. 



E-mail: kluOphysics . muni . cz 



1 Introduction 



One of the most important problem in present cosmology is the understanding of the 
origin of the late-time cosmic acceleration (the so-called Dark Energy (DE) epoch). 
Recently new interesting model DE model was proposed in [TJ |2]. This model 
consists of two scalar fields where one of scalars represents the Lagrange multiplier. 
The multiplier puts constraint on the second scalar field and as a result the theory 
contains singe degrees of freedom. It was shown that the energy of the system flows 
along time-like geodetic that is similar to the dust, however the theory contains non- 
zero energy. The behavior of this system suggests that it can be natural candidate 
for unification of Dark Energy and Dark Matter. The cosmological implications of 
these models were then analyzed in [H EJ E]. The role of Lagrange multipliers in 
the context of f(R) gravities was studied in [3]. Moreover, the Lagrange multipliers 
in the context of modified gravity may improve the ultraviolet properties of the 
covariant Hofava-Lifshitz gravity [9] leading to its renormalizability conjecture [3 [8]. 

As was shown in all these papers the presence of the Lagrange multipliers in 
the action has strong impact on the form of the resulting equations of motions. 
Then it is natural to ask the question how the presence of Lagrange multipliers 
modifies Hamiltonian structure of given theory. Moreover, we would like to see 
whether the Hamiltonian of these systems is again given as a linear combination 
of constraints and whether these constraints are the first class and their Poisson 
algebra respects the basic principles of geometro dynamics [TQ1 EU CE2] • It turns out 
that Hamiltonian structure of given theory is very interesting. We show that the 
presence of the first scalar field that plays the role of the Lagrange multiplier implies 
an existence of the second class constraints. Then after their solving we find the 
Hamiltonian equations of motions for the second scalar field that are autonomous in 
the sense that the time evolution of the scalar field does not depend on its conjugate 
momenta. Such systems were studied in the past especially in the context of the 't 
Hooft deterministic approach to quantum mechanics [131 El EE HB]- We also find 
that the resulting theory is a fully constrained system with the algebra of constraints 
that has the same form as in General Relativity. 

As the second example of the Lagrange multiplier modified theory we consider 
the gravity action introduced in [3]. This action is the Lagrange modification of 
F(R) gravity theories @. We show that the resulting Hamiltonian is given as a linear 
combination of constraints and has similar structure as the Hamiltonian of F(R) 
gravities (2D 122] I- 

However there is an important difference that follows the fact 
that the presence of the Lagrange multiplier implies that the original auxiliary fields 
become dynamical in Hamiltonian formulation. We further determine the Poisson 
brackets between constraints. We show that the algebra of these constraints takes 
exactly the same form as in [TOl [TT| IT2] . In other words we explicitly prove the 
consistency of Lagrange modified theories of gravity from the Hamiltonian point of 
view. 

2 For review, see [T7| [T9l [20] . 
3 For related works, see [24l [25] . 
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Let us summarize our results. We study the Lagrange multiplier modified the- 
ories with emphasis on their Hamiltonian formalism. We find that the resulting 
Hamiltonian is again given as a linear combination of the first class constraints. 
We show that the Poisson brackets of these constraints have the same form as in 
General Relativity. 

This paper is organized as follows. In the next section (J2j) we perform the 
Hamiltonian formulation of the General Relativity action together with the Lagrange 
multiplier modified scalar field action. We find corresponding Hamiltonian and 
diffeomorphism constraints and calculate their algebra. In section (J3]) we study 
the Lagrange multiplier modified action introduced in [3J. We again determine 
corresponding Hamiltonian. Then we calculate the Poisson brackets of the secondary 
constraints and we find that they take exactly the same form as in General Relativity. 

2 Lagrange Multiplier Modified Scalar Field Ac- 
tion 



In this section we develop the Hamiltonian formalism for Lagrange multiplier mod- 
ified scalar field action. We study the form of the action that was introduced in 


S = J S D+ ^xy^[^R(g) - ^^<R0 - V(cf>) - A[i^<R0 + [/(«/>)]] . 

(1) 

Let us explain our notation. We consider D + 1 dimensional manifold Ai with 
the coordinates x M , fi = 0, . . . , D and where x M = (t, x) , x = (a; 1 , . . . , x D ). We 
presume that this space-time is endowed with the metric g^ u (x p ) with signature 
(— , +,...,+). Suppose that M. can be foliated by a family of space-like surfaces 
Ei defined by t = x . Let gij,i,j = 1, . . . , D denotes the metric on Ej with inverse 
g l i so that gijg^ k = 8f. We introduce the future-pointing unit normal vector n M to 



the surface E 4 . In ADM variables we have n° = \/—g 00 , n % = —g° 1 / \/—g 00 - We also 
define the lapse function N = 1/ \J—g m and the shift function N l = —g 0t /g 00 - In 



terms of these variables we write the components of the metric g^ v as 
5-00 = — iV 2 + N i g lj N j , g 0i = N t , fa = g tj , 

9 N 2 ' 9 N2 , 9 9 N2 ■ 



(2) 



Then it is easy to see that 

V-detg = N^ffetg~ , g^d^A = ~(V n 0) 2 + g ij drfdjcj) . (3) 
Further, D + 1-dimensional curvature ( D+1 ^R can be written as 

(d+i )r = K i jK .. _ R 2 + R (D) + Zg^^/ZjjnVK) - -Z—diiy/ggVdjN) , (4) 

V 9 v 9^ 
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where = j^(d t gij — ViNj — VjiVj) and where V» is covariant derivative defined 
using the metric gij. Further, K = g^Kji. In what follows we ignore these boundary 
terms when we will presume appropriate boundary conditions. Let us consider the 
scalar field action. Using the notation introduced above we find the momentum 
conjugate to <p and A 

P<t> = Vg(u + A)V n , p x « . (5) 
Then the Hamiltonian for the scalar field takes the form 

H* = J d D *H^ , = N74 + N l nf , Hi = , 

(6) 

Finally we write the Hamiltonian for General Relativity part of the action 

H GR = J d D *H GR , H GR = NH GR + N l H GR , (7) 

where 

U GR = ^ 3 9^ kl -^ 2 -V~9R (D) , 
U GR = -2g lk V 3 Ti k \ 

(8) 

where n % i is momentum conjugate to g^ with non-trivial Poisson brackets 

{^(x), n k \y)} = + <5$)*(x - y) , (9) 

and where 7r = n^gji- Note that V, is a covariant derivative calculated with the 
metric g^ that also obeys Vigjk = 0. 

In summary, the total Hamiltonian is H = + H GR . The preservation of the 
primary constraints Pn ~ ,p l ~ implies the secondary ones 

u T = uf 1 + u% « o , Ui = u GR + nf^o . (io) 

It is useful to introduce the smeared form of these constraints 

T T (N) = T% R (N) + Tj.(iV) , 
T S (JV*) = T GR {N l ) + T£(7V*) , 

(11) 
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where 

T% R (N) = J d D *NH% R , T*(N) = J d D *NH\ 
Tf R (^) = J dPxrfHf* , T*(N l ) = J d^i^nf + N'pxdA) , 

(12) 

where we included the primary constraint p\ ~ into definition of T^iV*) in order to 
ensure the correct form of the Poisson bracket between the diffeomorphism generator 
T^(^) and the scalar field A. 

It is well known that the Poisson brackets between smeared form of the General 
Relativity constraints take the form [TUJ [TTJ H~2] 

{T% R (N),T% R (M)} = Tg R (g ij (NdjM — MdjN)) , 
{Tf fl (JV*),Tg B (M)} = T^i^diM) , 
{Tf^^.Tf^AT)} = Tjj^N'djM* - M'djN*) . 

(13) 

On the other hand we have to determine the Poisson brackets between constraints 
corresponding to the scalar field. First of all it is easy to see that 

{T^iV*), T*(Af )} = T (k s (N j d j M i - M j djN l ) . (14) 
On the other hand the Poisson bracket between Ts{N l ) and T^(M) is equal to 

'ts^n*), t£(m)} = J d D *(-N k d k n^ - d k N k n^) = 

d D xN k O k H* = T*(N k d k M) 



(15) 



using 



{T 5 (in^} = -N^gij-d^gkj-gikdje 
{T S (N*),^} = -N k d k ^g-^gd k N k . 



(16) 



Note that the presence of the term N i p\diX in the definition of Tg(iV*) was crucial 
for deriving of the correct form of the Poisson bracket f fl5|) . Finally we calculate the 
Poisson bracket between T^(iV), T^(M) and after some algebra we find the desired 
result 

{tJ(A0,tJ(M)} = T^iNdjM - MdjN)) . (17) 
It is also easy to show that 

T£*(iV),Tj(M)} + {t+(N),T% r (M)} = 

(18) 
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due to the fact that H T depends on g and not on their derivatives. If we combine 
these results we find that the Poisson brackets of the constraints TV(iV), Ts(iV*) 
has the desired form ffTS"]) . 

As the next step we analyze the stability of the primary constraint p\ ~ 0. The 
requirement of its stability implies the secondary constraint 

a tPA (x) = {p A (x), if} = 2 ^_ ( j +A)2 ^ - lv99 ij di(f>d j( j> -Vg~u = g x {*) « o . 

(19) 

We observe that 

fa(x), £ A (y)} = * + A)3 ^(x)^(x - y) . (20) 

In other words p\ and Q\ are the second class constraints. However there are addi- 
tional non-zero Poisson brackets. The first one is 

{&(x),&(y)} = -2 1 j^(x)#- y )- 



a. 



1 u ,r mx) 



*(x - y) . 

(21) 



It is also clear from the structure of the constraint Q\ that there is non-zero Poisson 
brackets between Q\ and Ti defined by ( TTUT) 



{£ A (x),?/}^0, (22) 

where T-L = NT-It + N l T-Li. Note that the explicit form of this Poisson bracket is not 
important for us. 

Using these results we can proceed to the study of the stability of the secondary 
constraints. Following the standard analysis of the constraint systems we introduce 
the total Hamiltonian as 

H T = H + J d D x(a Px + pg x ) , (23) 

where a, (3 are Lagrange multipliers and analyze the stability of the constraints 
Ti,G\,P\- Firstly we have 

W = {n(*),H} + J d D y(p{n(x),g x (y)} + Hn*)iPx(y)})^ 

d D y/3{H(x),£ A (y)}^0, (24) 



where we used the fact that the Poisson brackets between % and H weakly vanish. 
Then the requirement of stability of the constraint % ~ determines the value of 
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the Lagrange multiplier f3 to be equal to 0. On the other hand the time evolution 
of the constraint Q\ is given by the equation 

$&(x) = {Sa(x),#} + f d D ya(y) {£ A (x),p A (y)} » . 

(25) 

Due to the fact that {Qa,H} ^ and {G\,p\} ^ the equation above can be 
solved for a at least in principle. Then using these results it is easy to see that the 
constraint p\ ~ is preserved during the time evolution of the system. Further, p\ 
and Q\ are the second class constraints that can be solved for A and p\ so that the 
reduced phase space is spanned by (<7y, 7r lJ '), and the symplectic structure is 

given by the Dirac brackets between these variables. In order to find their form we 
introduce following notation for the Poisson brackets of the second class constraints 
Px,Gx 

A n (x,y) = Wx),p A (y)} = 0, A 12 (x,y) = {p A (x),£ A (y)}^0, 
A 21 (x,y) = {£ A (x),p A (y)}^0, A 22 (x, y) = {£ A (x), £ A (y)} ^ 

(26) 

and denote the inverse matrix as (A~ 1 ) AB (x, y). This matrix by definition obeys 
the equation 

J rf x A^(x,z)(A- 1 )^(z, y ) = ^(x-y) . (27) 
It can be shown that the matrix (A -1 ) has following structure 

=(:;)• (28) 

where * denotes non-zero elements. It is important for the calculation of the Dirac 
brackets that (A -1 ) 22 = 0. Explicitly, the Dirac bracket between and p^ takes the 
form 

{0(x),p (y)} D = {0(x),^(y)}- 

- J ^z^z'{0(x),$ A (z)}(A- 1 ) AB (z,z'){<I»B(z'),^(y)} = 

= {0(x),^( y )}-|^z^z'{0(x),^ A (z)}(A- 1 ) 22 (z,z'){^ A (z'),^(y)} = 
= {0(x),^(y)} , 

(29) 

where $a = (px, Qx) is the common notation for the second class constraints. 

We are now ready to completely eliminate the second class constraints The 
constraint Q\ = can be solved for u + A 

(u + A) = . n . (30) 
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Inserting this result into the Hamiltonian constraint ([6]) we find that it takes the 
form 

U% = p^^ + 2U+ - ^loU . (31) 

We observe that this Hamiltonian density is linear in momenta. Then the equation 
of motion for takes the form 

d t (t> = {0, H} = Ny/gvdiWrf + 2U + N%(j) (32) 

that shows that the time evolution of does not depend on p$. Such systems were 
extensively studied in the past in the context of 't Hooft's deterministic approach to 
quantum mechanics [131 US UHl US] and it is really interesting that the Hamiltonian 
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with similar structure arises in Lagrange modified multiplier theory. @ We complete 
our analysis by calculation of the Poisson bracket between "Hj, given in (131]) and the 
spatial diffeomorphism constraint Ts(iV*). Using 

{TsW^diWrf} = -N k d k {gvd^drf) (41) 

we easily find 

{^(iV 4 ), U(\ = -N*d#4 - diN*74 , (42) 

where 1-Lj, was given in (|3TT) . The analysis of the remaining Poisson brackets is the 
same as above with conclusion that the smeared form of the constraints obey the 



4 We review the basic facts considering such system, following [18]. Let us consider the Hamil- 
tonian system 

H =Pif i (q) + U(q) ,i = l,...,N. (33) 
From (13"3"|) we determine the equations of motion for q % 

d t q l = {q\H}=r(q) . (34) 

This equation for q % is autonomous, i.e., it is decoupled from the conjugate momenta pi. Further 
it is impossible to perform the Legendre transformation to the Lagrangian since ify — d p.g p . = 0. 
However it is possible to find the Lagrangian that gives the equation of motion (|3"4"]l when we 
introduce the auxiliary fields A^ and write the Lagrangian as 

L = X t (q'-f(q))-U(q). (35) 

No we show that from ([33)) we can derive the Hamiltonian (|3"3")l . The momenta conjugate to A,; 
and q % take the form 

A = wr ' * ! = S= A - (36) 

so that we have two sets of primary constraints 

$\=p\n0, $f=pf-A l «0. (37) 
The extended Hamiltonian that follows from (1351) takes the form 



H E = H + tj^i + lo^I , H = Xif + U(q) . (38) 



Then we study the stability of the constraints , $ 
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d t $\ = {$\,H E } = -f i + uji = 

dp 
dq 



(39) 



From these equations we can in principle determine the Lagrange multipliers. In other words the 
constraints $^ > are the second class that should strongly vanish. The solving of these constraints 
we find the Hamiltonian 

H=p q l P + U{q) (40) 



that coincides with the Hamiltonian (|33|) . Further, it can be easily shown that the Dirac brackets 
between q l and pi coincide with their Poisson brackets. However the problem with the Hamiltonian 
([3"3")) is that is not bounded from below which is due to the absence of a leading kinetic term 
quadratic in the momenta (pt) 2 . 
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algebra of constraints given in ffl3|) . In other words we show that the Lagrangian 
multiplier modified scalar action together with General Relativity action obeys the 
basis rules of geometrodynamics. 



3 Hamilton Analysis of F(R) Theories with La- 
grange Multipliers 

It turns out that the Lagrange multiplier modified _F(.R)-gravity possesses many 
interesting properties. For example, the reconstruction programme can be more 
easily performed in Lagrange multiplier modified gravity [3]. In usual F(i?)-gravity, 
we need to solve the complicated differential equation to realize the reconstruction 
program, for recent review, see [23]. It was demonstrated in [3] that the presence 
of constraint significantly simplifies the reconstruction scenario. It was also shown 
there that the presence of Lagrange multiplier implies that it is necessary to include 
the second F(R) function into action. 

The action introduced in [3] takes the form 



S = I d D+1 xy/-g 



F^P+VR) - A f^d^Rg^d^R + F 2 (^R)^j 



(43) 

Introducing two auxiliary fields A, B we can rewrite the action (14*31) into the form 

S = J d D+1 xy/^g~ F X {A) - A (^d^Ag^d v A + F 2 (A)^j + B^ D+1 ^R - A) . (44) 

It is easy to see that integration of A, B from (T44|) leads to (j4"3l . Our goal is to find 
the Hamiltonian from ( 14*4")) implementing D + 1 formalism. In fact using (jlj) it is 
easy to see that the action (jHJ) takes the form 

S = J d D xdt^N (F 1 (A) - X(-W n AW n A + g ij diAd 3 A + F 2 {A)) - BA) + 
+ J d D xdt^NB(K ij Q ljkl K kl + R (D) — A) — 
- -2 J d D xdt(y/g(d t B — N i d i B)K + 2y/gdiBg ij djN) , 

(45) 



where we performed integration by parts and ignored boundary terms. From (1431) 
we easily find momenta conjugate to canonical variables g^, N,N iy A and B 

7r« = JgBQ^Ku - JgV n Bg« , ,p N ^0, J w , 
p B = -2^K , p A = 2V\V n A , p x w . 

(46) 
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Note that the Lagrange multiplier implies that A is a dynamical field which is 
different from standard F(R) theory of gravity where A remains auxiliary field. 
Then after some effort we derive the Hamiltonian density in the form 



where 



14 1 ^ij- „ JA 1 2 ^Pb 

ti T = J 9ikguK 'FTTn + 

+ ^=p{D-l)p l B -JgBlt D ) + 2d i [ y /gg i id j B] 

+ -QjyPA + VdBA - y/g[Fi(A) - XisTdiAdjA + F 2 (A))} 



and where 



= -^p\ + ^BA-^[F 1 {A)-\{g^d i Ad j A + F 2 {A))\. 



(47) 



(4£ 



-Hi = PA diA + p B diB + pxdA - 2g ik V j n jk . (49) 
For further purposes we split %t into two parts as Ht = + where 

T ^ B y*y* yfgBD ^gD 

+ I ^{D-l)p 2 B - y /gBR^+2d i [ y /gg ij d j B) 



(50) 



The theory possesses 2 + D primary constraints 

TX N « ,7Ti « ,7T A « . (51) 

The preservation of the primary constraints itn and 7Tj imply the secondary con- 
straints %t ~ , ~ while the preservation of rr\ « leads to the secondary 
constraint 

^ A = " ^(^diAdjA + F 2 (A)) » (52) 

We see that it takes the same form as the secondary constraint ( fT9l) . Clearly p A 
together with are the second class constraints. Properties of these constraints 
were analyzed in previous section and results derived there can be used in this 
section as well. 

On the other hand the form of the Hamiltonian constraint Hj, R is new and we 
have to check that this constraint is preserved during the time evolution of the 
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system. In other words we have to calculate the Poisson brackets of the smeared 
form of these constraints 

T$ R (N) = J d D xiV(x)H^(x) , T G S R (N' 1 ) = J d D xiV i (x)Hf i? (x) . (53) 

Let us now outline the strategy of the calculations of these Poisson brackets. In 
the process of their calculations several delta functions occur. However it turns out 
that the non-zero contributions give terms that contain derivatives of these delta 
functions. Such expressions arise for example from following Poisson bracket 

K . (x ),(^) (y)} = -^W). (54) 



The right side of this equation can be calculated using the formulas 
5RW = -( R (D)yi5 gi . + - g ij V k V k 6 9ji , 5g = gg ij 5 9ij 



(55) 



Now we are ready to perform these calculations. It turns out that following non-zero 
Poisson brackets contribute to the final result 



J d D xN^-^g ikgjl i^,J d D yM^gR^ 



j d D yNy/gR( D \J d D ^M^-=7r ij g ik9jl ^ 



= 2 J d D x(NViVjM - MViVjN)7i ij - 2 j d D ^(NV i V i M - MV { V l N) + 
+ A J d D x.Tc ij (NViM — MViN)ir ij -^\7 jB — 4 J d D x7c(NV \M - MV 'iN)^V*B 



(56) 



= -— J d D ^ii(NV i V i M - MViVW) + 2 j d D ^{NV l V i M - MV^N) - 
- — / d D xn{NViM - MVtN)—^- + 4 / rf D XTr(iVV i M - MV { N) 



V { B 



B 

(57) 



5 In [22] similar analysis has been performed in the context of non-projectable version of Hofava- 
Lifshitz F(R) gravity. 
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and 

{/ dDxN ^I dD y^ R{D)BM ] + {/ d D yVgR {D) BN,J AM-^|} = 

^ _ D " / dD *P BB ( NViViM - M ViVW) - 
y d D xp B (NViMM - MViN^B , 

(58) 



2(1 -£) 
D 



d D ^^^g lk g^ k \ J dy2Md i [ y /gg ij d j B]^ + 
+ dy2Nd i [^gg ij d j B]J d D ^^^g lk g^ kl J = 



2 y rf D x7r-^(A^V i M-MV i AT)^ V jJ B-4 y cPx-^TW^M - MViN^VjB 

(59) 



and 



y dD *^=BD* 2 >f d D yM2d m [^g mn d n B] 
d D xN2d m [^g mn d n B], J d D x-^n 2 > = 



2(2 -D) f jD n 



D 



J d D x^(NViM - MViN)V l B 



(60) 



and 



jg D *PB,j d D y2Md m [^gg mn d n B]^ + 
d D y2Nd rn [^g mn d n B], J d D x-^np B 
(2 ~ D) y rf D x(iVV m M - MV m N)p B V m B + 



+ |j y d D X7r(NV m V m M — MV m V m iV) , 



(61) 
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d°x- NB 



. ^-1)^,2 J d D yMd m [^gg mn d n B)} + 
+ |2 y d D xNd m [^g mn d n B], J d D y^^(D-l)pl 



D - 1 
D 



J d D *p B B{NV m V m M - MV m V m A) . 



(62) 



Collecting all these terms together we obtain that almost all contributions cancel 
and the final result takes the form 

{T% R (M),T% R (N)} = Tl R ((NVjM - MVjN)g ji ) . (63) 

In other words the Poisson bracket of the smeared form of the Hamiltonian con- 
straints fl50l) has the same form as in General Relativity and hence it is with agree- 
ment with basic principles of geometrodynamics. Alternatively, it has the form that 
is expected for fully diffeomorphism invariant theory. Note also that the Poisson 
bracket between smeared form of the diffeomorphism and Hamiltonian constraint 
takes the standard form that follows from the fact that Hamiltonian is manifestly 
invariant under spatial diffeomorphism. Then it is clear that the diffeomorphism 
and Hamiltonian constraints are preserved during the time evolution of the system. 

Now it is straightforward to finish the analysis of the Poisson brackets of the 
constraints of the Lagrange multiplier modified gravity. Since the Poisson brackets 
of the constrains corresponding to the gravity part of the action are the same as in 
General Relativity and since the scalar part of the constraints has exactly the same 
form as in previous section we immediately find that the Poisson brackets of the 
Lagrange multiplier modified F(R) gravity take the form 

{TY(A0,T T (Af)} = TstfUNdjM-MdjN)) , 
{T 5 (iV*),T r (M)} = T T (N i d i M) , 
{T 5 (A^ i ),T 5 (M i )} = T s (N j d j M i - M'djN*) . 

(64) 

where 

n T = U% R + v.? , Hi = ~g a W k Ti lk + PA diA , (65) 
where H^ R is given in (HB1) . Note that Ti^ is equal to 



•Ht = VA\j g ij diAdjA + F 2 (A) + ^BA - y/gF x {A) (66) 

after solving the second class constraint Q\ given in ( 152]) with respect to A 

A = PA . (67) 

2 y /g y /F 2 (A)+g i W i (l)d j <f> 
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